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ARTIFACT REMOVAL FROM AN ELECTRIC SIGNAL 
Field of the Invention 

The present invention relates generally to a method and apparatus for 
5 suppressing or removing artifacts from an electric signal. The term "artifact" 
here refers to interference which has a certain structure and is therefore not 
noise-like random interference. The invention is intended especially for 
medical technology to eliminate artifact components from an 
electrocardiogram (EGG) signal in connection with magnetic resonance 
10 imaging (MRI), but the method can be applied to any systems in which an 
artifact component similar to the one caused by an MRI device to an EGG 
signal is superimposed on a useful signal. Artifact components are typically 
caused by sources which are external to the apparatus receiving the useful 
signal. 

15 

Background of the Invention 

As is known, MRI (Magnetic Resonance Imaging) is a technique used 
for getting images of a body without the use of X-rays. In MRI, the patient is 
placed within a rapidly changing magnetic field created by the MRI device. 

20 Radio waves are simultaneously transmitted to the patient, and images of the 
interior of the body are obtained by means of a computer analysis of the radio 
waves received from the patient. 

During the imaging, the condition of the patient is constantly 
monitored. This includes the monitoring of the EGG signal of the patient, for 

25 example. One drawback related to the monitoring is that the intense magnetic 
fields of the imaging device induce interference with the EGG signal of the 
patient. This Interference, which is termed here an artifact, is short-termed and 
strong, and appears very frequently. 

FIG. la and lb show a typical example of the narrow and strong 

30 interference peaks caused to the EGG signal by MRI measurement. FIG. la 
illustrates a sequence of the narrow interference pulses superimposed on the 
EGG signal, while FIG, lb illustrates the shape and length of one interference 
pulse. The artifact is usually rather similar for each MR image acquired and 
consists of a few short peaks. The shape of an interference peak is roughly 

35 box-like, with a relatively constant period between the steepest rise and fall 
edges. 

The prior art methods intended for suppressing the above-described 



artifacts are based on either limiting the rise and fall speed of the ECG signal 
or freezing the signal during the interference periods. In the following these 
prior art methods are discussed briefly. 

A natural way of detecting the MRI gradient peaks is to monitor the 
5 rise and/or fall speed, i.e. the slew rate, of the ECG signal. A solution based on 
a slew rate limiter is disclosed in U.S. Patent 4,991,580. This patent presents a 
method whereby the ECG signal is supplied to the input of a slew rate limiter 
(SRL) circuit, which limits the slew rate of the ECG to a pre-selected maximum 
value slightly greater than a typical maximum slew rate in an ECG signal. 

10 Another approach for eliminating MRI induced artifacts superimposed 

on an ECG signal is presented in U.S. Patent 5.038,785. This patent describes 
a cardiac monitor in connection with an MRI device, the monitor to include a 
comparator which compares each wave form received from the ECG 
electrodes with the properties of a cardiac signal, such as the slope. When the 

15 comparator detects interference in the ECG signal, it gates a track and hold 
circuit. The track and hold circuit passes the ECG signal, except when gated 
by the comparator. When gated by the comparator, the track and hold circuit 
continues to supply the same output voltage it had at the beginning of the 
gating period, i.e. it freezes the ECG signal for the period of interference. 

20 U.S. Patent 4,243,045 presents a method for the suppression of 

narrow interference pulse peaks and of interference hum in a useful signal. 
The method is intended particularly for suppressing the pacemaker pulses and 
the power supply hum superimposed on an ECG signal. The narrow pulses 
* caused by the pacemaker are blanked out by keeping the signal constant for 

25 the blanking-out period, and the power supply hum is eliminated by an 
interference hum filter. The filtered-out interference hum component is 
opposingly superimposed on the original signal to compensate for the 
interference hum. 

The prior art methods based on freezing the signal entail the 
30 fundamental drawback of not being able to take into account signal changes 
occurring during the freezing period. The methods based on limitation of the 
slew rate, in turn, involve the fundamental drawback that some of the 
interference remains in the useful signal, as the methods limit only the size of 
the interference. For these reasons, the prior art methods may leave errors of 
35 substantial magnitudes affecting the useful signal. 

It is the objective of the invention to improve the accuracy of the 
methods described above and to bring about a new and efficient solution for 



3 

suppressing or eliminating artifacts from a useful signal, particularly MRI- 
induced artifacts from an ECG signal. 

Summary of the Invention 

5 The objective of the present invention is to bring about a solution for 

eliminating as accurately as possible narrow interference peaks from a useful 
signal, such as MRI-induced artifacts from an ECG signal. 

According to the invention, the useful signal is first sampled at a high 
sampling frequency. The samples are then supplied to an artifact detection 

10 process during which the degraded samples are Identified, i.e. it is detected 
which of the samples are degraded by the artifact and which are essentially 
free of artifact and thus good samples. The selected high quality samples of 
the useful signal are then utilized in a subsequent downsampling process, 
which is in this context temried an artifact-aware downsampling, since the 

15 preceding detection process indicates the samples containing an artifact 
component. More particularly, in the downsampling process, samples 
considered as free of artifact are used to generate a new sequence of samples 
having a sampling frequency substantially lower than the original sampling 
frequency. By utilizing the artifact-free samples only, the downsampling 

20 process thus constructs a new sample sequence with a reduced sampling 
rate. Consequently, the method of the invention includes a controlled 
downsampling of the original signal, the downsampling being controlled by the 
artifact-content of the high rate sample sequence. 

Thus one aspect of the invention is providing a method for 

25 suppressing interference in an electric signal, the method comprising the steps 
of 

- sampling the electric signal at a first sampling frequency, whereby a 
first sequence of samples is obtained, 

- selecting some of the samples in the first sequence of samples on 
30 the basis of predetermined criteria, and 

- downsampling said first sequence using the selected samples, 
whereby a second sequence of samples is obtained, the second sequence 
forming a digital presentation of the electric signal in which the interference is 
suppressed. 

35 The method of the invention improves the accuracy of artifact 

elimination, since it does not limit or freeze the signal values. Instead, the 
signal is sampled at a high rate, the artifact-content of the samples is detected. 
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and the detected artifact content controls the downsampling of the signal so 
that only high quality sannples are utilized. 

Another aspect of the invention is that of providing an apparatus for 
suppressing interference in an electric signal, the apparatus connprising 
5 - sampling means for sampling the electric signal at a first sampling 

frequency, whereby a first sequence of samples is obtained, 

- selection means for selecting some of the samples in the first 
sequence on the basis of predetermined criteria, and 

- downsampling means, operatively connected to the selection 
10 means, for downsampling the first sequence by using the selected samples, 

whereby a second sequence of samples is obtained, the second sequence 
forming a digital presentation of the electric signal in which the interference is 
suppressed. 

In a preferred embodiment of the invention, the samples of the first 
15 sequence are compared with a reference signal generated from the original 
useful signal, whereby the selecting step is performed based on the comparing 
step. 

In a further preferred embodiment of the invention, information about 
the artifact content of the samples is generated as a result of the comparison 

20 of the signal with the reference. The information generated is supplied, 
together with the samples, to the downsampling process, where the signal is 
downsampled using this information. The information consists preferably of a 
sequence of flags, each flag being attached to the relevant sample and 
indicating where the sample is artifact-free or not. 

25 In another preferred embodiment of the invention, the power supply 

hum, possibly appearing at 50 or 60 Hz and superimposed on the ECG signal, 
is eliminated from the ECG signal before the downsampling process. This is 
implemented by first calculating the amount of power supply hum in the ECG 
signal by using the artifact-free samples only. The calculated value of the hum 

30 at each sample is then deducted from the corresponding sample, whereby 
updated sample values are obtained for the high rate sample sequence. In this 
way the power supply hum appearing in the signal can be efficiently removed. 
Thus one advantage of the method of the invention is that it allows 
simultaneous and efficient suppression of both MRI induced artifacts and 

35 power supply hum. 

Other features and advantages of the invention will become apparent 
by reference to the following detailed description and accompanying drawings. 
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Brief Description of the Drawings 

In the following, the invention and its preferred embodiments are 
described more closely with reference to the examples shown in FIG. 2 to 3c 
5 in the appended drawings, wherein: 

FIG. la illustrates the artifact superimposed on an ECG signal in 
connection with magnetic resonance imaging (MRI), 

FIG. 1b illustrates one artifact pulse in FIG. 1a, 

FIG. 2 illustrates the apparatus according to the invention. 
1 0 FIG. 3a illustrates one embodiment of the method of the invention, 

FIG. 3b illustrates a second embodiment of the method of the 
invention, and 

FIG. 3c illustrates a third embodiment of the method of the invention. 

1 5 Detailed Description of the Invention 

FIG. 2 is a schematic presentation of the apparatus according to the 
present invention. The signal received from the ECG sensor(s) is supplied to 
an amplifier stage 210, which amplifies the signal before it is sampled and 
converted into digitized format in an A/D converter 220. A high sampling rate is 

20 used in the A/D converter to sample the ECG signal. A typical sampling rate 
could be of the order of 10 kHz, for example, which yields about 10 samples 
during one interference peak (cf. FIG.Ib). In any case, the initial sampling rate 
used in the present invention is higher than the sampling rate needed for the 
original signal, the latter being typically 100 to 1000 Hz in the case of an ECG 

25 signal. 

As a result of the sampling, a high rate sample sequence MRS is 
obtained, which is buffered and processed in a digital signal processing unit 
230. The processing unit performs the removal of the artifacts and outputs the 
corrected ECG signal, i.e. the original ECG signal from which the artifacts have 

30 been removed. 

FIG. 3a illustrates one embodiment of the method performed in the 
digital signal processing unit 230. A separate reference signal is first formed 
from the high rate sample sequence HRS (step 310). This reference signal is 
preferably generated from the sample sequence itself, through the use of a 

35 slew rate limiter, for example. In other words, the reference signal is preferably 
generated by copying the input signal and limiting the slew rate in the copy. 
The reference signal serves as a reference for identifying the artifact in the 
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original input signal. 

Instead of being generated from the input signal, the reference signal 
may also be an external signal, i.e. a signal which is not derived from the ECG 
signal input to the signal processing unit. The MRI device, for example, can 
5 provide an extefnal reference signal whose high amplitude levels identify the 
samples containing an artifact component. In view of the invention, the actual 
method of generating the reference signal is not essential. 

The samples of the input signal containing an artifact component are 
then identified (step 320). Although these samples can be identified by 

10 comparing an external reference signal, such as the above-mentioned signal 
provided by the MRI device, to a predetermined threshold, it is assumed here 
that the said samples are identified by comparing the input signal, i.e. the 
original ECG signal, with the reference signal generated from the original ECG 
signal by limiting the slew rate. The comparison is performed by comparing 

15 each sample of the input signal to the corresponding sample of the reference 
signal. When a sample of the input signal meets predefined criteria, it is 
considered as an artifact-free sample. Different criteria may be used in 
determining when a sample is artifact-free. In this example, when the input 
signal and the reference signal are within the noise level of each other, the 

20 sample is considered artifact-free and thus valid. 

Based on the comparison process at step 320, some of the high rate 
samples are selected and the selected samples are forwarded for a 
downsampling process. In the embodiment of FIG. 3a, the selection is 
implemented as follows. First, on the basis of the results of the comparison 

25 process, information is generated about the validity of each sample of the input 
signal. This information preferably consists of a sequence of flags, one flag 
being attached to each sample of the input signal, to indicate whether the 
corresponding sample is valid (i.e. free from an artifact). Second, the samples 
and the associated flags are forwarded for the downsampling process, during 

30 which the digital input signal is downsampled utilizing the information about the 
validity of the samples, i.e. the flag sequence (step 340). As the information 
about the artifact content of the samples, i.e. the flag sequence, is forwarded 
for the downsampling process, the process is here called artifact-aware 
downsampling. Based on the information, only artifact-free samples are 

35 selected, and the downsampling is performed based on these samples only. 

It is to be noted here that it is preferred not to discard the samples 
affected by artifacts before entering the downsampling process, but instead all 
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the samples are supplied for said process, together with the flag sequence, so 
that the downsampling process receives the information about the locations in 
the time domain of the good and bad samples. Another alternative is to discard 
the samples which are identified as having been affected by artifacts and to 
5 supply only the good samples for the downsampling process. However, this 
alternative requires that the significant moments corresponding to the good 
samples are forwarded for the downsampling process. 

The actual downsampling can be performed in many ways, which are 
known in the art. For integer downsampling ratios, the most common approach 

10 is to use an anti-aliasing low-pass filter and then to select every n*^ sample. For 
non-integer downsampling ratios, the signal is usually first upsampled by 
padding the signal with zeros, the signal is then low-pass filtered, and then the 
desired samples are selected (every n^ sample). To keep the individual 
sample rate changes small, the whole process is often divided up into several 

15 stages with smaller sampling rate ratios. 

Another approach used for downsampling is interpolation: the signal 
at the correct point in time is calculated by modeling the neighboring values. 
Methods typically used are to select the value of the nearest neighbor or to 
model the signal using a polynomial function or spline. 

20 As a result of the downsampling, a new sample sequence with a 

lower sample rate is obtained. The sample rate of this sequence can be 300 
Hz, for example. The new sequence forms the digital presentation of the ECG 
signal from which the artifact has been removed. 

FIG. 3b illustrates a preferred embodiment of the method performed 

25 in the digital signal processing unit. In addition to the functions described in 
connection with FIG. 3a, this embodiment includes the elimination of power 
supply hum, i.e. elimination of the interference at 50 or 60 Hz caused by a 
power supply. In the artifact-aware downsampling described above, each 
output sample is based on input samples having very small artifact levels. 

30 However, since the good samples are taken from varying parts of the time 
window, the output may have a disturbance resembling temporal jitter in the 
sampling. This will deform the signal slightly, whereby the signal can interfere 
with a notch filter intended to remove the 50/60 Hz power supply hum. 
Therefore, better results can be achieved by using the solution of FIG. 3b, 

35 where the power supply hum is removed before downsampling (step 330). The 
elimination can be implemented by estimating the amplitude and phase of the 
power supply hum and removing the estimated hum from the input signal, i.e. 
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subtracting the estimated hum from the original signal. The estimation is 
implemented by using only the samples considered as good ones when 
calculating the amplitude and phase of the power supply hum. Thus in the 
process of removing power supply hum, the information about the artifact 
5 content of the samples is utilized, similarly as in the downsampling process. 
The power hum estimate can be calculated by means of a band-pass filter, 
which is thus an artifact-aware filter as it uses the good samples only. 

When the downsampling process comprises a FIR filter, the artifact- 
aware downsampling can be obtained by calculating the weighted averages of 

10 only the good samples with their normal weights and by correcting the gain 
based on the weights of the bad samples. For example, if the FIR filter is a 
boxcar filter, the artifact-aware downsampling process, i.e. the filter, calculates 
the average of the corresponding good samples, said average to form the 
value of one low rate sample. 

15 When the downsampling process includes interpolation, the standard 

interpolation means can be applied by taking into account the good samples 
and the corresponding moments in the time domain. 

In a further preferred embodiment of the invention, the above- 
mentioned artifact-aware downsampling is combined with artifact-aware low- 

20 pass filtering. The filter can be replaced by an artifact-aware FIR filter using the 
number of good samples as the flag. This kind of implementation further 
improves the accuracy of the apparatus, and it may also be more tolerant to 
long periods of noise. 

Low-pass filtering may be a separate step subsequent to 

25 downsampling, or it may be divided into two stages, of which the first one is 
formed by the above-mentioned low-pass filtering in the downsampling 
process. The second stage may then occur after the downsampling, whereby 
the filtering in the second stage is applied to the lower rate sample sequence. 
The division of low-pass filtering into two stages is advantageous in terms of 

30 calculating power; the processing time is shorter if the filtering is divided into 
two stages of which the first one provides rough filtering and the second one is 
a more accurate stage. The use of low-pass filtering is illustrated in FIG. 3c. 
The low-pass filtering stage 350 of receiving the low rate sample sequence 
and the corresponding low rate flag sequence can thus form the only low-pass 

35 filtering stage of the apparatus or it can complement the low-pass filtering 
performed at stage 340. 

By using more than one bit in the flag, an indication may be given of 
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the magnitude of the artifact contained in the sample. In this way the artifact 
components can be classified into several categories depending on their 
magnitude. As the quality of a low rate sample calculated depends on the 
number of high rate samples used to produce it, the flag sequence supplied for 
5 stage 350 can be numeric, and the FIR filter can be implemented by using all 
the samples but having different weights for each sample. 

The artifact-aware FIR filter may be implemented according to 
equation (1), for example: 

i 

10 where x is the input signal, b is the filter coefficierit vector, and f is the 

flag sequence vector. A binary flag sequence is represented with zeros for bad 
samples and with ones for good samples. 

The artifact-aware band-pass filter producing the estimated power 
supply hum sequence at step 330 may be implemented according to equation 

15 (2), for example: 

yp = A(FArFx (2) 

where the elements of the nx2 matrix A are Ai,i = sin(27ti/n) and Ai.2 = 
cos(27ri/n), F is a diagonal matrix having the flag sequence as the diagonal 
elements, {-f is the Moore-Penrose pseudoinverse, and x is the input signal. 

20 Although the invention was described above with reference to the 

examples shown in the appended drawings, it is obvious that the invention is 
not limited to these, but may be modified by those skilled in the art without 
departing from the scope and spirit of the Invention. Even if the method was 
described with reference to the MRI-induced artifact superimposed on the 

25 ECG signal, the method can be applied to any system in which an artifact 
component with like properties is superimposed on a useful signal. 



